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Linear  Time-Variant  Space-Variant  Filters,  Wave  Propagation  In  a  Random 
Medium,  The  WKB  Approximation,  Time-Space  Transformaticas,  Output  Electrical 
Signals,  Coherence  Function 
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Wave  propagation  In  a  random.  Inhomogeneous  ocean  Is  treated  as  transmission 
thru  a  linear,  time-variant,  space-variant,  random  comaunlcatlon  channel.  A 
consistent  notation  (vla-a-vls  ad  hoc),  fundamental  Input-output  relations, 
and  various  time-space  transformations  for  both  deterministic  and  random 
linear,  time-variant,  space-variant,  filters  are  established.  Using  the 
method  of  separation  of  variables  and  the  W.K.B.  approximation,  a  time- 
invariant,  space-variant ,  random  transfer  function  of  the  ocean  volume  Is 
derived.  The  ocean  volume  Is  characterized  by  a  random  Index  of  refraction 
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•  function  of  depth.  The  index  of  refrectlon  la  decomposed  Into  a 
deterministic  component  and  s  zero  mean  random  component.  In  addition,  two 
example calculetiona  are  made.  The  first  example  involves  the  derivation  of 
the  equations  for  the  random,  output  electrical  signals  at  each  element  in  a 
planar  array  of  complex  weighted  point  sources  In  terms  of  the 
frequency  spectrum  of  the  transmitted  electrical  signal,  the  transmit  and 
n^nye,  and  the  transfer  function  of  the  ocean  medium.  The  second 
example  Involves  the  derivation  of  the  coherence  function,  l.e.,  the  auto¬ 
correlation  function  of  the  transfer  function. 
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ABSTRACT 


Wave  propagation  in  a  random,  inhomogeneous  ocean  is 
treated  as  transmission  thru  a  linear,  time-variant,  space- 
variant,  random  communication  chemnel.  A  consistent  notation 
(vis-a-vis  ad  hoc) ,  fundeuaental  input-output  relations,  and 
various  time-space  transformations  for  both  deterministic  and 
random  linear,  time-variant,  space-variant, filters  are  esta¬ 
blished.  Using  the  method  of  separation  of  variables  and  the 
W.K.B.  approximation,  a  time-invariant,  space-variant,  random 
transfer  function  of  the  ocean  volume  is  derived.  The  ocean 
volume  is  characterized  by  a  random  index  of  refraction  which 
is  a  function  of  depth.  The  index  of  refraction  is  decomposed 
into  a  deterministic  component  and  a  zero  mean  random  compo¬ 
nent.  In  addition,  two  example  calculations  are  made.  The 
first  example  involves  the  derivation  of  the  equationfeVor 
the  random,  output  electrical  signals  at  each  element  in'a 
receive  planar  array  of  complex  weighted  point  sources  in 
terms  of  the  frequency  spectrum  of  the  transmitted  electrical 
signal,  the  transmit  and  receive  arrays,  and  the  transfer 
function  of  the  ocean  medium.  The  second  example  involves 
the  derivation  of  the  coherence  function,  i.e.,  the  auto¬ 
correlation  function  of  the  transfer  function. 
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I.  INTRODUCTION 

Since  the  wave  equation  for  small  amplitude  acoustic  sig¬ 
nals  is  linear,  we  can  represent  the  ocean  medium  as  a  linear, 
time-variant,  space-variant,  random  filter  (system  or  communi¬ 
cation  channel)  in  general.  With  this  interpretation  in  mind, 
refer  to  Fig.  1  which  illustrates  a  basic  bistatic  communica¬ 
tion  channel  geometry,  and  Fig.  2,  which  is  a  mathematical 
block  diagrcun  representation  of  Fig.  1.  With  respect  to  Fig.  1, 
both  the  transmit  and  receive  apertures  (arrays)  are,  in  gen¬ 
eral,  volume  apertures  (arrays)  and  in  motion.  Before  pro¬ 
ceeding  further,  a  word  of  caution  concerning  Fig.  2.  Note, 
for  example,  that  ^  ^  ^  general. 

The  equations  required  for  describing  the  filter's  input-output 
relationships  and  for  coupling  the  transmitted  and  received 
electrical  signals  to  the  medium  via  the  transmit  and  receive 
apertures  are  developed  in  Sections  II  and  III,  respectively. 

Let  us  now  describe  the  notation  used  in  Fig.  2.  The  posi¬ 
tion  vectors  r  and  r  refer  to  spatial  coordinates  (x  ,y  ,z_) 
and  (x,y,2),  respectively,  and  t  refers  to  time  in  sec.  The 
parameters  f  and  n  are  frequencies  in  HZ.  where  f  represents 
input  or  transmitted  frequencies  while  i  represents  output  or 
received  frequencies.  Note,  that  if  n  ^  f,  Doppler  spread  is 
implied. 

The  quantities  a,  £,  and  are  vectors  whose  components 
are  spatial  frequencies  with  units  of  cycles/m.  Since  spatial 
frequencies  are  related  to  both  direction  cosines  and  wavelength. 
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and  hence,  wavenumber  components,  they  represent  directions 
of  wave  propagation.  The  vector  y_  represents  input  or  trans¬ 
mitted  spatial  frequencies  into  the  medium  as  in  X^^Cf,^), 
while  £  represents  output  or  received  spatial  frequencies 
from  the  medium  as  in  Y„(n,S).  Note,  that  if  3  ^  ,  angular 

spread  (scatter)  is  implied. 

The  remaining  expressions  found  in  Fig.  2  are  further 
described  in  the  following  list: 

x(t,r)  -  input  electrical  signal  to  transmit  electro¬ 
acoustic  transducer  applied  at  time  t  and 
spatial  location  r  of  transducer. 

X(f,^)  -  frequency  (f)  and  angular  (o_)  spectrum  of  input 
electrical  signal. 

A.j,(f,r)-  complex  frequency  response  at  spatial  location 
r  of  transmit  transducer.  Also  referred  to  as 
the  complex  transmit  aperture. 

D^(f,a)-  transmit  far-field  directivity  function  or  beam 
pattern. 

Xj^(t,r)-  input  acoustic  signal  to  the  medium  applied  at 
time  t  and  spatial  location  r.  Also,  output 
acoustic  signal  from  transmit  electro-acoustic 
transducer. 

Xjj(f,v)-  frequency  (f)  and  angular  (^)  spectrum  of  input 
acoustic  signal. 

hjj(T,r^;  t,r)  -  time-variant,  space-variant  impulse  re¬ 
sponse  of  the  ocean  medium.  It  represents  the 
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response  of  the  laedlum  at  time  t  and  spatial 
location  r  due  to  the  application  of  an  unit 
impulse  at  time  (t-T}sec.,  or  t  sec.  ago,  at 
a  distance  lr-r^|  m.  away  (see  Fig.  1). 

Hj^(f t,r)  -  time-variant,  space-variant  transfer 
function  of  the  ocean  medium. 

yj^(t,r)  -  output  acoustic  signal  from  the  medium  at  time 
t  and  spatial  location  r.  Also,  input  acoustic 
signal  to  receive  electro-acoustic  transducer. 

-  frequency  (n)  and  angular  (S)  spectrum  of  out- 
put  acoustic  signal. 

Aj^(n,£)  -  complex  frequency  response  at  spatial  location 
r  of  receive  transducer.  Also  referred  to  as 
the  complex  receive  aperture. 

-  receive  far-field  directivity  function  or  beam 
pattern. 

y(t,r)  -  output  electrical  signal  from  receive  electro¬ 
acoustic  transducer  at  time  t  and  spatial  loca¬ 
tion  r  of  transducer . 

~  frequency  (n)  and  angular  (y)  spectrum  of  output 
electrical  signal. 

As  was  mentioned  previously,  we  can  represent  the  ocean 
medium  as  a  linear,  time-variant,  space-variant,  random  filter. 
The  term  "time-variant"  implies  motion  ^unongst  targets,  the 
ocean  surface,  discrete  point  scatterers,  and  the  transmit  and 
receive  apertures  (arrays) .  Discrete  point  scatterers  in  the 
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ocean  nay  include,  for  example,  gas  bubbles,  fish,  and  other 
particulate  matter.  The  time-variant  property  results  in 
both  Doppler  spread  zmd  spread  in  round-trip  time  delay  values. 
If  the  filter  is  time- invariant,  then  no  motion  is  implied. 

As  a  result,  there  will  be  no  Doppler  spread  and  no  spread  in 
round-trip  time  delay. 

The  term  ” space-variemt*  implies  that  the  sound  speed 
profile  (index  of  refraction)  of  the  ocean  is  a  function  of 
position.  The  space-variant  property  results  in  scatter  or 
angular  spread  due  to  refraction.  If  the  filter  is  space- 
invariant,  then  an  isospeed  medium  is  implied.  As  a  result, 
there  will  be  no  refraction,  and  hence,  no  scatter  or  angular 
spread  since  the  sound  rays  will  be  travelling  in  straight 
lines. 

In  addition,  since  any  motion  and/or  the  index  of  refrac¬ 
tion  can  be  decomposed  into  a  sum  of  deterministic  (average) 
and  random  (fluctuating)  components,  these  random  components 
can  be  accounted  for  via  a  random  filter  representation  vis-a- 
vis  a  deterministic  filter  representation. 

By  using  a  systems  theory  approach,  surface,  volume,  and/or 
bottom  reverberation  returns  can  be  modelled  as  the  outputs 
from  linear  filters.  In  addition,  target  returns  can  also  be 
modelled  as  filter  outputs.  Furthermore,  different  transmit 
signals  and  transmit  and  receive  directivity  functions  can 
easily  be  coupled  to  various  models  (i.e.,  transfer  functions) 
of  the  random,  inhomogeneous  ocean  medium  in  a  straightforward 


and  logical  fashion  in  order  to  study  their  effects  on  target 
detection  or  parameter  estimation  using  various  space-time 
signal  processing  algorithms. 

The  approach  of  treating  the  ocean  as  an  isospeed  medium, 
and  hence,  as  a  linear,  time-variant,  random  communication 
channel  is  well  established  [1-19].  This  linear,  time-varying, 
random  system  theory  approach  has  also  been  applied  to  target 
scattering  problems  in  radar  astronomy  [20]  and  to  communica¬ 
tion  channels  in  general  [21-23].  However,  with  respect  to 
target  models,  past  research  efforts  have  been  devoted  mainly 
to  the  slowly  fluctuating  point  target  problem  [24-31].  Ef¬ 
forts  to  treat  more  complicated  target  models  were  made  by 
Kooij  [32],  Moose  [10],  and  Ziomeh  and  Sibul  [19,33].  Kooij 
[32]  and  Moose  [10]  both  modelled  the  target  as  a  Inear,  time- 
invariant,  deterministic  filter  while  Ziomek  and  Sibul  [19,33] 
modelled  the  target  as  a  linear,  time-varying ,  random  filter. 

In  addition,  Ziomek  [34]  has  shown  that  the  form  of  the  gen¬ 
eralized  ambiguity  function  can  be  derived  by  treating  the 
scattered  acoustic  pressure  field  from  a  point  target  (in 
relative  motion  with  respect  to  a  bistatic  transmit/receive 
array  geometry)  as  the  output  of  a  linear,  time-varying,  random 
filter . 

Some  work  has  been  done  in  treating  the  ocean  medium  as  a 
linear,  time-variant,  space-variant ,  random  filter  by  Laval 
[9,35]  and  Laval  and  Labasque  [36].  However,  the  notation 
used  to  incorporate  the  space-variant  property  is  ad  hoc,  i.e., 
spatial  variables  are  simply  included  in  the  arguments  of  the 
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impulse  response  and  transfer  functions,  for  exeunple,  rather 
than  having  evolved  from  a  systematic  and  consistent  notation 
based  upon  linear,  time-varying,  space-varying  system  theory. 

In  addition,  Laval  and  Labasque  136]  assume  functional  forms 
for  the  ocean  transfer  function  instead  of  deriving  them. 
Middleton  [37,38]  also  studied  underwater  acoustic  propagation 
in  a  random,  inhomogeneous  ocean,  but  did  not  concern  himself 
directly  with  the  derivation  of  random,  time-variant,  space- 
variant  ocean  transfer  functions.  He  described  the  propagation 
phenomena  using  space-time  operators. 

In  this  paper  we  will  study  underwater  acoustic  propagation 
in  a  random,  inhomogeneous  ocean  by  treating  the  ocean  medium 
as  a  linear,  time-variant,  space-variant,  random  filter.  Sec¬ 
tion  II  is  devoted  to  a  discussion  of  the  fundamentals  of 
linear,  time-variant,  space-variant  filters  and  is  based  upon 
the  generalization  of  the  results  contained  in  Ziomek  [39] . 

A  consistent  notation  is  developed  in  a  systematic  manner.  Va¬ 
rious  system  functions  are  introduced  and  important  input-output 
relations  and  multidimensional  (time-space)  Fourier  transform 
pairs  are  derived  for  both  deterministic  and  random  filters. 

To  the  best  of  the  author's  )tnowledge,  the  expressions  presented 
in  Section  II  have  not  appeared  previously  in  the  literature. 

The  equations  necessary  to  couple  the  medium's  transfer 
function  to  the  far-field  becun  patterns  of  the  transmit  and 
receive  apertures  (arrays)  and  to  the  frequency  spectrum  of 
the  transmitted  signal  are  discussed  in  Section  III. 
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Section  4.1  is  devotad  to  the  main  problem  of  the  paper 
which  is  the  derivation  of  a  random  ocean  transfer  function 
incorporating  the  W.K.B.  approximation.  In  the  process  of 
the  derivation,  the  index  of  refraction  is  decomposed  into 
deterministic  and  random  components  and  it  is  assumed  that 
the  medium  is  "weakly  scattering".  The  transfer  function 
derivation  was  motivated  by  the  work  of  Clarke  [40]  on  the 
application  of  the  W.K.B.  approximation. 

Section  4.2  is  devoted  to  an  example  calculation  of  the 
equations  for  the  random,  output  electrical  signals  appearing 
at  each  element  in  a  receive  planar  array  of  complex  weighted 
point  sources  in  terms  of  the  frequency  spectrum  of  the  trans¬ 
mitted  electrical  signal,  the  transmit  and  receive  arrays,  and 
the  random  ocean  medium  transfer  function. 

Finally,  in  Section  4.3,  the  autocorrelation  function  of 
the  transfer  function,  which  is  also  known  as  the  coherence 
function,  is  calculated. 
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II.  FUNDAMENTALS  OF  LINEAR  TIME-VARIANT  SPACE-VARIANT  FILTERS 
2. 1  Deteraiinistic  Filters 

2.1.1  Impulse  response  and  transfer  functions 
A  linear,  time-variant,  space-variant,  filter  is 
depicted  in  Fig.  3,  where  it  is  characterized  by  its  corres¬ 

ponding  time-varying,  space-varying  impulse  response  h(T,r^; 
t,r) .  The  function  h(T,r^;  t,r)  describes  the  response  of  the 
filter  at  time  t  and  spatial  location  r  =  (x,y,z)  due  to  the 
application  of  an  unit  impulse  at  time  (t-t) ,  or  t  seconds  ago, 
and  at  a  distance  |£-£q|  meters  away  where  r^  =  ^’‘o'^o'^o^’ 

Note  that 

h(T,r  ;  t,r)  i  h(t,r;  t-t,  r-r  ) .  (2.1-1) 

— O  —  —  —o 

The  relationship  between  the  input  signal  x(t,r) 
and  the  output  signal  y(t,r)  is  given  by 


y(t,r) 


r-r  )h(T,r  ;  t,r)dTdr 
—  — o  — o  —  — o 


(2.1-2) 


where  it  should  be  noted  that  both  the  input  and  output  signals 
are  functions  of  time  and  space. 


Exeunole  2.1-1 

Different  forms  of  Eq.  (2.1-2)  can  be  obtained  by  ma)cing 
the  following  simplifying  assumptions: 
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(1)  if  the  linear  filter  h  ia  time- invariant  and  space-inva- 
riant,  then 


h(T,r^;  t,r)  =  h(t  -  [t-x] ,  r  -  [r-r^] ) 

=.  h(t,r^)  (2.1-3) 

and  as  a  result,  Eq.  (2.1-2)  reduces  to 


y (t,r) 


r-r  ) h (x , r  ) dxdr 
—  — o  — o  — o 

(2.1-4) 


which  is  a  multidimensional  convolution  integral  as  would 
be  expected. 

(2)  if  h  is  time-invariant  but  space-variant,  then 


h(x,r^;  t,r)  «  h(t  -  [t-x],  r;  £-£q) 


»  h(T,r;  r-r^) 

-  h(x,r  ;  r)  (2.1-5) 

— o  — 


and  as  a  result,  Eq.  (2* 1-2}  reduces  to 


y(t,r) 


m 

//, 


;e<t-x. 


r-r  ) h ( T , 
—  — o 


r  ;  r)dTdr  . 
-o  -  -o 


(2.1-6) 
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(3)  if  h  is  time-variant  but  space-invariant,  then 

h(T,r^;  t.r)  =  h(t,  r  -  Ir-r^] ;  t-r) 

=  h(t,r^;  t-T) 

=  h(T,r^;  t)  (2.1-7) 

and  as  a  result,  Eq.  (2.1-2)  reduces  to 


y(t. 


OB  00 

-// 


x(t-T,  r-r^)h(T,r  ;  t) dtdr. . 
—  -o  -o  — o 


(2.1-8) 


Note  that  if  Eq.  (2.1-2)  is  rewritten  as 


y (t,r) 


a,  r-i)h(a,£;  t,r)dod£ 


(2.1-9) 


and  if  the  input  is  an  unit  impulse  applied  at  time  a  =  (t-x) 
seconds  at  a  distance  |^|  =  Ir-r^|  meters  away,  i.e.,  if 

x(a,5_)  *  J(a  -  [t-rj,  £  -  [£-r^]  )  .  (2.1-10) 

then  substituting  Eq.  (2.1-10)  into  Eq.  (2.1-9)  yields 
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h(T,r  ;  t,r) 
— o  — 


F"^F“^{HCf,v.;  t,r)l 


C2.1-14) 


or 


h(T,r^;  t,r) 


t,r)exj.  > 


j  2it  ft )  exp  (- j  2it^*  r^)  df d^. 


(2.1-15) 


The  choice  of  a  plus  (+>  sign  in  the  exponent  of 
exp (  + j 2iTy^- r^)  appearing  in  Eq.  (2.1-13),  which  corresponds  to 
the  forward  spatial  Fourier  transform  w.r.t.  r^,  was  not 
arbitrary.  This  choice  of  sign  convention  is  meant  to  be  con¬ 
sistent  with  that  of  the  spatial  Fourier  transform  relation¬ 
ship  between  a  complex  aperture  function  and  its  directivity 
function  (beam  pattern)  as  is  developed  later  in  Section  3. 
Besides,  the  integrand  term 

exp(-j2irfT)exp(  +  j2Ttv*r^) 

appearing  in  Eq.  (2.1-13)  has  the  nice  physical  interpretation 
of  being  a  time-harmonic  plane  wave  travelling  in  the  direction 
of  increasing  r^  =  |r^l  with  the  sign  convention  as  given. 

Example  2.1-2 


If  the  linear  filter  h  is  time- invariant  and  space-in¬ 
variant,  then  its  corresponding  transfer  function  can  be  ob- 
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talned  by  substituting  Eq.  C2-1-3)  into  Eq.  t2.1-13).  Doing 
so  yields  the  following  result; 


H(f ,y) 


m  ea 

If 


hCr  ,r^)exp<-j2iifT)exp(+ j2ity^«r^)dTdr^ 


(2.1-16) 


Similarly,  from  Eq.  (2.1-15), 


h(T,ro) 


Ok  oa 

J  j 


)  exp(+j2itfT)exp(-j2nv^*r^)dfdv^. 


(2.1-17) 


Calculating  the  transfer  function  according  to 
Eq.  (2.1-13)  assumes  that  one  Icnows  the  impulse  response  func¬ 
tion.  However,  even  if  the  impulse  response  is  not  known, 
the  transfer  function  can  still  be  obtained  by  using  a  time- 
harmonic  plane  wave  as  an  input  signal.  This  can  be  shown  by 
representing  a  linear,  time-varying,  space-varying,  filter  by 
the  linear  operator  L(>)  which  operates  on  input  signals  that 
are  functions  of  both  time  and  space.  The  output  of  the 
filter  y(t,r)  can  then  be  expressed  as 


y (t,r)  »  L(x(t,r) ] . 


(2.1-18) 


» 
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If  we  let 


x(t,r)  “  exp(+j2itft)exp(-j2itv'r)  ,  (2.1-19) 

which  is  in  the  form  of  a  time-harmonic  plane  wave  travelling 
in  the  direction  of  increasing  r  =  |rj,  then  from  Eqs.  (2.1-2) 
and  (2.1-18)  we  can  write 


L  [exp(+j2)tft)exp(-j2iiv-r)] 


■»  CD 

II 


exp[  +  j2iif  (t-r)]exp[-j2itv*  (£-£_)  1 


h(T,r  ;  t,r)dTdr  , 


(2.1-20) 


and  by  using  the  definition  of  the  transfer  function  given  by 
Eq.  (2.1-13),  Eq.  (2.1-20)  can  be  expressed  as 

L[exp  (+ j  2  ff  t)  exp  (- j  2Trv  •£ )]  =  H(f,jj;  t,  r)  exp  (  + j  2Tf  t)  exp  (- j  2ii\^*  r) 

(2.1-21) 

where  L  [exp  (  + j  2itf t)  exp  (-j2n^*r)  ]  is  the  response  of  the  filter 
to  exp(+j2irft)exp(-j2irv"£),  and  H(f,^;  t,£)  is  the  time-varying, 
space-varying  transfer  function  of  the  filter  evaluated  at 
f  and 

Equation  (2.1-21)  is  a  fundamental  result  that 


will  be  used  in  the  derivation  of  an  ocean  transfer  function 
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2.1.2  Spreading  and  bi-frecruencv  functions 
Two  additional  filter  functions  will  now  be  intro¬ 
duced;  namely,  the  spreading  function  and  the  bi-freguencv 
function.  The  spreading  function  will  be  discussed  first. 

The  spreading  function  S(T,r^;  ♦,£)  is  defined  as 


i(t,r^;  ♦,<.)  A  p^F^{h(T,r^;  t.r)  }  (2.1-26) 


or 


m  m 

l(t,r^;  ♦,«.)  ^  J  j h(T,r^;  t,r)exp(-j2»*t)exp(+j2iiic.r 


)dtdr 
(2.1-27) 


where  corresponds  to  the  rate  of  change  of  the  filter's  im¬ 
pulse  response  in  HZ.  and  ^  is  a  vector  whose  components  are 
spatial  frequencies  which  correspond  to  the  rate  of  change  of 
the  impulse  response  in  cycles/m.  Similarly, 

h(T,r  ;  t,r)  *  FT^Fr^{S(T,r  ;  ♦,<)}  (2.1-28) 


or 


h(x,r^;  t,r) 


♦  ,£)  exp(+j2it»t)exp(-j2iric_*£)d^dic_. 


(2.1-29) 
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An  alternate  representation  of  the  output  y(t,r) 
can  be  obtained  by  substituting  Bq.  v.2.1-29)  into  Eg.  (2.1-2). 
Doing  so  yields 


y(t,r) 


m  cm  m  m 

fin 


X  (t-t,  j  -’'♦t)  exp  (- j  2Tf^-r ) 


S(T,r  ;  ♦,  pc)d4>d<dTdr 


(2.1-30) 


where  the  integrand  term 


S(T,r  ;  ♦,<)x(t-T,  r-r  )exp(+j2ii4it)exp(-j2if<T) 

is  a  time,  space,  and  "frequency”  (frequency  in  HZ.  and  spatial 
frequencies  in  cycles/m.)  shifted  replica  of  the  input  signal 
x(t,r)  weighted  by  the  spreading  function  S(T,r^; 

Therefore,  for  a  given  input  (source)  location  as  specified  by 
r^,  the  spreading  function  determines  the  amount  of  spread  in 
round-trip  time  delay  1(360.),  frequency  ♦(HZ.),  and  spatial 
frequencies  j<  (cycles/m.)  that  an  input  signal  will  undergo  as 
it  passes  through  a  linear,  time-varying,  space-varying,  channel. 
Equations  (2.1-30)  and  (2.1-2)  both  indicate  that  the  output 
y(t,r)  at  some  receiver  location  r  is  dependent  upon  the  loca¬ 
tion  r^  of  the  input  (source)  via  the  spreading  function 
S(T,r^;  ♦,£)  w.r.t.  Eq.  (2.1-30),  and  via  the  impulse  response 
h(T,r^;  t,r)  w.r.t.  Eq.  (2.1-2).  For  example,  the  existence  and 
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extent  of  SOFAR  transmission  depends  on  both  the  depths  of  the 
source  and  receiver  [41J. 

The  last  filter  function  to  be  discussed  is  the 
bi-frequency  function.  The  bi-frequency  function  B(f,y^;  if,^) 
is  defined  as 


3(f,u;  ♦,<)  ^  P,.Fj.{H(f,v;  t,r)  } 


(2.1-31) 


or 


OB 

l(f,y.:  ♦»<)  ^  //  H  (f  t,  r)  exp(- j2iitt)  exp  (+ j2iT£*  r)  dtdr . 


(2.1-32) 


Similarly, 


H(f,v;  t,r)  »  F^^F~^{B(f,v;  ♦,<)} 


(2.1-33) 


H(f,v:  t,r)  = 


eg  SB 

// 


B  (f  ,y^;  ♦ ,  <.)  exp(+ j2:t4)t)  exp  (- j 2nK_*  r )  d*dK_. 

(2.1-34) 


The  term  bi-frequency  function  was  originally  used  to  denote 
the  appearance  of  the  two  frequency  variables  f  and  ♦  in  linear, 
time-varying  systems  theory.  We  can  generalize  this  notion 
by  using  the  term  bi-frequency  function  to  denote  the  appearance 
of  the  input  "frequency"  pair  f  and  y^,  and  the  filter  variation 
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"frequency"  pair  ^  and  ^  in  linear,  time-varying,  space-vary¬ 
ing,  systems  theory. 

Just  as  the  spreading  function  gives  an  indication 

of  how  rapidly  the  impulse  response  changes  with  time  and 

space,  the  bi-frequency  function  gives  an  indication  of  how 

rapidly  the  transfer  function  changes  with  time  and  space. 

The  bi-frequency  function  can  also  be  obtained  by 

taking  the  Fourier  transform  of  the  spreading  function  w.r.t. 

T  and  r  ,  i . e . , 

— o 


B(f,vj  4i,<)  =  F  F  {S(T,r;  <>,«:)} 

_  _  X  _ 


(2.1-35) 


or 


B ( f ,  v ;  i  ,<) 


0 ,  <_)  exp  (- j  2i  f  :  )  exp  (*32-’_^-r^)d:dr 

(2.1-36) 


I 


and,  similarly, 


S(T,r^;  t,<)  =  F"^F"^{B(f , 5,^) 


(2.1-37) 


s  (  r ,  £0 ;  ii ,  <_) 


5  ,£)  exp  (  + j  2nf  I )  exp  (- j  2:iy^-r  )dfdy^  . 


(2.1-38) 


I 


In  addition,  the  impulse  response  and  bi-frequency 
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2.1.3  Output  frequency  and  angular  spectrum 
The  output  frequency  (n)  and  angular  (^)  spectrum 
Y(n,B)  is  defined  as 


Y(n,3)  =  F^F^{y(t,r) > 


(2.1-43) 


Y(n,a) 


•  OB 

// 


y  ( t,  r)  exp  (- j2'»nt)  exp  (+ j  r)  dtdr  (2.1-44) 


where  n  corresponds  to  output  frequencies  in  HZ.  and  £  is  a 
vector  whose  components  are  output  spatial  frequencies  with 
units  of  cycles/raeter .  Substituting  Eq.  (2.1-25)  into  Eq. 
(2.1-44)  yields 


Y(n,3) 


OS 


woe  woe 


t,r) exp[- j2T (i-f ) t]  • 


exp [+j 2tt  (£-^)  'rjdtdrdfd^ 


(2.1-45) 


where,  from  Eq.  (2.1-32),  it  can  be  seen  that  the  inner  multi¬ 
dimensional  integral  is  equal  to  3(f,v;  n-f,  £-y)  so  that 


Y(n,S) 


woe 


n-f  ,  3-\; )  dfdy 


(2.1-46) 


I 
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Which  is  in  the  form  of  a  multidimensional  "frequency"  domain 
convolution  integral.  Note  that  ♦  =  n-f  or  n  =  t+t,  and  that 
K  5  8-^  or  8.  5  Thus,  the  output  frequencies  n  in  HZ.  are 

equal  to  the  sum  of  the  input  frequencies  f  and  the  variations 
in  frequency,  ♦,  due  to  the  time-varying  property  of  the  filter. 
Similarly,  the  output  spatial  frequencies  i  (directions  of 
wave  propagation)  in  cycles/m.  are  equal  to  the  sum  of  the  in¬ 
put  spatial  frequencies  ^  (directions  of  propagation  of  trans¬ 
mitted  plane  waves)  and  the  variations  in  spatial  frequency 
(directions  of  wave  propagation) ,  jc,  due  to  the  space-varying 
property  of  the  filter.  Equation  (2.1-46)  demonstrates  that  a 
linear,  time-variant,  space-variant,  filter  will  spread  the 
input  frequency  and  angular  spectrum  X(f,v)  in  both  frequencies 
in  HZ.  and  spatial  frequencies  in  cycles/ra. 

Example  2.1-3 

If  the  linear  filter  h  is  time- invariant  and  space-in¬ 
variant,  then  H(f,^;  t,r)=  H(f,^)  [see  Eq.  (2.1-16)]  and,  as  a 
result,  Eq. 

B(f,^;  ♦,jc) 

or 

B(f,v;  ^,K)  =  H(f,v)6(4i)6(i)  (2.1-47) 


(2.1-32)  reduces  to 


QB 

H(f,v)  j 


/' 


exp(-j2':T^t)dt  /  exp(  + j27Tjc-£)dr 


I 


L 
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since 

F^{1}  =  4{$1 

and 

F^{1)  =•  S{<)  . 

Substituting  Eq.  (2.1-47)  into  Eq.  (2.1-46)  yields 


V(n.3) 


X(f  ,v)H(f , v) 5 (n-f) 5 (8-v)dfdv 


which  simplifies  to 

y(n,3)  =  .x(n,3)H(n,3) , 

and  by  replacing  n  with  f  and  ^  with  we  finally  obtain 

Y(f,v)  =  X(f,v)H(f,v)  (2.1-48) 

which  is  the  output  frequency  and  angular  spectrum  from  a  linear, 
time- invariant,  space- invariant ,  filter.  Note  that  the  output 
frequency  f  and  the  output  spatial  frequency  vector  ^  are  iden¬ 
tical  with  those  of  the  input  XCf.y^).  Hence,  as  would  be  ex¬ 
pected,  there  is  no  "frequency"  spreading.  This  would  corres¬ 
pond  to  the  physical  situation  of  transmitting  a  signal  via  a 
transmit  aperture  (array)  to  a  receive  aperture  (array)  when 
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the  platforms  containing  the  apertures  (.arrays!  are  not  in 
motion,  and  the  intervening  ocean  medium  has  a  constant  speed 
of  sound  (index  of  refraction)  and  no  discrete  point  scatterers 
in  motion.  No  motion  implies  no  Doppler  (frequency)  spread 
and  a  constant  speed  of  sound  implies  that  the  sound  rays 
travel  in  straight  lines,  i.e.,  there  is  no  angular  spread 
(scatter)  or  change  in  the  direction  of  propagation  of  the 
transmitted  sound  rays. 


***** 
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2 . 2  Random  Filters 

2.2.1  Filter  autocorrelation  functions 
In  Section  2.1,  four  system  functions  were  intro¬ 
duced  which  are  used  to  characterize  linear,  time-variant, 
space-variant,  filters.  However,  if  the  filter  is  random, 
then  each  of  these  system  functions  must  be  considered  as 
random  functions.  As  a  result;,  we  must  work  with  the  respec¬ 
tive  system  autocorrelation  functions  which  are  defined  as 
follows : 

Rh • T ' .r  1^1;  t,t',r,r')  *  E(h(T,r  ;  t,r) h* (t ' .r* ;  t',r')} 
n  o  “  o  ““  ■“  o  "" 

(2.2-1) 

Rjj{f,f',;i,v';  t,t',r,r')  »  E(H  (f ,  v;  t,£)  H*  (f  ' ,  v  ’ ;  t',r’)} 

(2.2-2) 


4  E{S(T,r  ; 

-o 


♦,<)S* (t • 

—  — o 


(2.2-3) 


and 

Rg(f,f  '  v' ;  ♦.o'.jc.k')  ^  E(B(f,v;  »  ,  <  )  B»  (  f  '  ,  v  '  ;  $',<’)} 

(2.2-4) 

where  £(•}  is  the  expectation  or  ensemble  average  operator 
and  the  asterisk  denotes  complex  conjugate. 

If  we  use  the  sign  convention  that  forward  trans¬ 
forms  w.r.t.  T,t,r^,  and  r'  are  defined  with  a  minus  sign  in 
the  exponent  of  the  complex  exponential  (inverse  transforms 
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w.r.t.  and  are  defined  with  a  plus  sign)  and  for¬ 

ward  transforms  w.r.t.  t*  t',r^,  and  r  are  defined  with  a 
plus  sign  (inverse  transforms  w.r.t.  f','f>'.^»  and  £  are  defined 
with  a  minus  sign) ,  then  it  can  be  shown  that  the  four  system 
autocorrelation  functions  are  related  by  the  following  Fourier 
transform  pairs: 


Rjj(f,f,^,v';  t.f.  r,r‘) 

and 

t.f  ,£,£•) 


^^•^r  t,f.£,£')} 

— o  — o 

(2.2-5) 


F*^F"l-F~^F"l-{Rg(f  ,f  •  ,  v.u*  ;  t.f  ,£,£’)  ); 


(2.2-6) 


.Rg(f  ,f  '  ,  V,  v' ;  .  {Rjj  (f.  f :  t.f, £,£')} 

(2.2-7) 

and 


(2.2-8) 


^  ' 'So'^'  ♦.♦'.f.'f.'J  =  ^t^t'^r^r  ’  ^^h^*  '  ^  t.f  ,£,£')} 

(2.2-9) 

and 


t,f,£,£’)  =  F"^F'|F'^F^J-{Rs(T,T',£^,f  ; 


-0-0 
(2. 2-10) 


Rgff 


,f’,£,v';  >  f,<,^>)  =  F  F..F  F  {R  (t,t',£^,£';  s , j ■,<,<■) ) 

T.ro£g3  -0-0  - 


(2.2-11) 
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(2.2-12) 


and  finally. 


Rg(f,f ',  v,w';  ♦,♦•,£,£')  =•  F^F^.Fr  F^'FtFfFi-Fr' 

t.fr.r')}  (2.2-13) 


=  F;^F;"F;V^F;iF;^F;Vy{R3(f.f  ,v,v', 


(2,2-14) 


The  interdependence  which  exists  amongst  the  four 
filter  autocorrelation  functions  is  illustrated  in  Fig.  5. 


dence  amongst  the  four  filter  autocorr* 
that  characterize  linear,  time-varying 
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2.2.2  uncorrelated  spreading  and  the  scattering 
function 

In  this  section  we  will  examine  the  consequences 
of  assuming  that  the  spreading  function  S(r.r^;  is  un- 

correlated  with  S(f,  r^;  for  ili  values  of  f  ^  t, 

r'  r  ,  and 

— O  “O 

The  assumption  of  uncorrelated  spreading  is  mathe¬ 
matically  equivalent  to  stating  that  the  autocovar Lance  of 
and  S<T',r^;  ♦',£')  is  aero  for  ^  values  of 

t'  ^  T,  ♦'  ^  ♦'  ^  '  i •  ® • ' 


Cg(T,T'  .rg.t^;^.^' 


S 

r'^r,  ♦’?<>»  and  < ' 

^  r  f 

(2.2-15) 


where  Cg  is  the  autocovariance 


function  and  ti£) 


E{S(t,L3J  ♦,£)}.  If  it  is  assumed 
then  Eq.  (2.2-15)  is  equivalent  to 


assumed  that  Wg(Tr£Q'’  ♦»£.)  * 


R  (t  T',r  /r ' :  j  <  •  K ' )  *  Rglt^r^;  ♦»£.)  5  ' )  * 


5 (r  -r* ) ^(♦-♦’ ) 4 (£“£' )  (2.2-16) 

-o  -o 


where 


Rg(T,Lj;»,£)  *E(lS(T,rQ;  ♦#£,)  I  ) 


(2.2-17) 
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is  called  the  scattering  function  and  is  equal  to  the  mean 
squared  value  of  the  spreading  function.  The  scattering  func¬ 
tion  can  be  thought  of  as  an  average  power  density  function 
which  determines  the  average  amount  of  spread  that  an  input 
signal's  power  will  undergo  as  a  function  of  round-trip  time 
delay  t  >  frequency  ^ ,  and  spatial  frequencies  £  for  a  given 
input  (source)  location  Note  that  the  scattering  function 

is  a  real,  positive  valued  function. 

Equation  (2.2-16)  is  the  result  of  the  assumption 
that  the  spreading  function  is  zero  mean.  However,  if  the 
spreading  function  is  non-zero  mean,  it  is  convenient  to  do 
the  analysis  with  the  centered  process 

S^(T,r  ;  ♦,<)  =>  S(T,r  ;  ♦,<)  -  Uc(T,r  ; 

C  — o  “  — O  —  a  — o  — 

The  random  process  S^(T,r^;  has  zero  mean,  and  as  a  re¬ 

sult,  its  autocovariance  function  is  equal  to 


C_  (T,T',r  ,£■;  =  E{S^(T,r  ;  ♦ , « )  S*  (*  '  ,£•  r  ')'',<’)} 

— O  —o  —  —  C  — O  c  — o  — 


Cs(T,t' ,r^,r^, 


i.e.,  the  autocorrelation  function  of  the  ’ero  mean  centered 
process  is  equal  to  the  autocovariance  funu.icn  of  the  original 
non- zero  mean  spreading  function. 
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Let  us  next  examine  the  effect  that  the  assumption 
of  uncorrelated  spreading  has  on  the  remaining  three  system 
autocorrelation  functions.  As  Fig.  5  indicates,  can 
be  obtained  from  Rg  by  performing  the  following  transformations: 


Rjj(f,f.v,.'  ;t,f, £.£•). 


(2.2-18) 


or 


Rjj  ?  t,t 


«•  c» 

'.£.£•)  =  I  ••• 


Rq  ( T  ,  T  * ,  ,  r^  t 


exp(+j2*((|)t-(>’t‘)]exp[-j2ff  (<.•£-<.’  •£')  I* 
exp(-j2ii  (f  T-f '  T  ' )  ]  exp[+j2ir  N'r^-y,’  -r^) )  • 
d(td<(i 'dicdic 'dtdt 'dr  dr  *  .  (2.2-19) 


If  Eq.  (2.2-16)  is  substituted  into  Eq.  (2.2-19),  then 


Rjj(f,f',y,y';  t,t*,£,£')  =  Rjj(Af,av:  At,A£)  (2.2-20) 


where 


R„(Af  , At,Ar) 


«  e»  e» 
//// 


♦  ,K)exp(+j2n$At)  • 


exp(- j  2:r£*  a£)  exp  (- j  2iiaf  t  )  • 
exp  (-•- j2TtA;^*r  )d^d<dTdr  , 


« 

* 

I 


(2.2-21) 
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Af  =  f-f,  4^  =  ^~}L' '  ~  t-t',  and  Ar  =  r-r '  . 

It  can  be  seen  from  Eq.  (2.2-20)  that  when  uncorrelated  spread 
ing  is  assumed,  the  autocorrelation  function  Rjj  becomes  a  func 
tion  of  the  differences  Af,  A^,  At,  and  Ar  only.  This  implies 
that  the  random  process  t,r)  is  wide-sense  stationary 

in  frequency,  spatial  frequencies,  time,  and  space.  An  addi¬ 
tional  requirement  for  H(f,^;  t,r)  to  be  wide-sense  stationary 
is  that  the  mean  value 

yjj(f,^;  t,r)  =  E{H(f,y;  t,r)  } 

is  a  constant.  Since  the  four  filter  functions  are  related 
by  linear  transformations  (see  Fig.  4)  ,  and  since  it  was 
assumed  that  Uc(T,r^;  ♦,<)  =  0,  then  u.  (T,r_;  t,r)  =  0, 
yjj(f,;^!  t,r)  =  0,  and  PgCf,^;  »,<)  =0  which  are  constant, 
zero  mean  values.  Therefore,  the  condition  of  uncorrelated 
spreading  in  round-trip  time  delay  input  (source)  location 
r^,  frequency  and  spatial  frequencies  <_  is  equivalent  to  a 
condition  of  wide-sense  stationarity  in  frequency  Af,  spatial 
frequencies  A^,  time  At,  and  output  (receiver)  location  Ar, 
respectively.  If  uncorrelated  spreading  in  x ,  r^,  » ,  and  ic 
occur  together,  then  we  have  a  wide-sense  stationary  uncorre¬ 
lated  spreading  (WSSUS)  communication  channel. 

Consider  the  autocorrelation  function  next. 

If  Eq.  (2.2-16)  is  substituted  into  Eq.  (2.2-10)  and  the  in¬ 
dicated  transformations  are  performed,  then 
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VT.x'.r  ,£■;  t,t*,r,r’)  =  at,  Ar)  5  (t-t  ' )  5  (r  -r" ) 

— o  —  n  —o  —  — o  — o 


(2.2-22) 


where 


Rj^(T,r^;  at,Ar) 


Rg{t,r^;  4>,£)exp(+j2iii)>At)  • 
exp  (-j2ir£^*Ar)djid£. 


(2.2-23) 


Equation  (2.2-22)  indicates  that  the  random  process  h(T,r^;  t,r) 

is  wide-sense  stationary  in  tine  and  space  because  of  the  it 

and  Ar  dependence  and  since  t,r)  =»  0.  Equation  (2.2-22) 

also  indicates  that  h(T,r^;  t,r)  is  uncorrelated  for  all  values 

of  T '  ■:  and  cl  r. 

— o  — o 

Finally,  if  Eq.  (2.2-lS)  is  substituted  into 
Eq.  (2.2-11)  and  t.he  indicated  transformations  are  performed, 
then 


Rg  (  f  ,  f  '  ,  y '  ;  ^,i}',y,y')  =  Rg  ( if ,  Ay;  5  ,  y)  5  ( $-•>  '  )  5  (y-y '  ) 

(2.2-24) 


where 


Rg ( if , Ay;  9 , < ) 


ji,y)exp(-j2xAfT)- 


exo(+j2xAu>r  )dTdr  . 
‘  —  — o  — o 


(2.2-25) 


Equation  (2.2-24)  indicates  that  the  random  process  3(f,y:  s.y) 
is  wide-sense  stationary  in  frequency  and  spatial  frequencies 


I 


I 
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because  of  the  4f  and  dependence  and  since  ♦,£)  =  0. 

Equation  12.2-24)  also  indicates  that  4,0  is  uncorre¬ 

lated  for  all  values  of  4*  /  4  and  jc'  ft  jc. 

Therefore,  in  summary,  under  the  asstnnption  of 
uncorrelated  spreading  the  four  filter  autocorrelation  func¬ 
tions  reduce  as  follows: 

4, 4', <»>«')  =  4.  If )  5  (t-t  • )  4  (r  -rl  )• 

4(4-4*) «(<-<’)  (2.2-16) 

Rjj  (f ,  f  *  V  ' ;  t,  t*  ,r,r ' )  =  Rjj(Af,A;y;  At,Ar)  (2.2-20) 

Rj^(t,t' ,r^,r^;  t,t',r,r')  =■  Rj^(T,r^;  At,  Ar)  4  (t-t  ' )  4  (r^-r^) 

(2.2-22) 

and 

Rg  (f ,  f  '  ,  V  ;  4, 4', £,<_')  =*  Rg  ( Af ,  A^;  4,  O  4  (  4-*  '  )  4  (£-£' ) 

(2.2-24) 

where  Rg(T,r^r  4,jf)  is  the  scattering  function. 


i 


I 
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2.2.3  The  scattering  function  and  its  Fourier 
transforms 

By  inspecting  Eq.  (2.2-23),  we  can  write  that 


I.,-!, 


R.  (x.r  ;  At, dr)  =  F  ;  >,<)> 

n  — o  —  o  <  o  — o  — 


(2.2-26) 


and 


or 


(2.2-27) 


A  ga 


I  h‘" 


Rs(’'£o  '  “  I  I  (~j2^>it)  • 

exp  (+ j  2jtK  •  Ar)  ditdAr .  (2.2-28) 


Also,  by  inspecting  Eq.  (2.2-25),  we  can  write 


that 


and 


Rg (Af , Av 


,  Ay  ;  ♦,<_)  =  ^ 


(2.2-29) 


Rc  ^  ^  f  3^ 

S  '-0 


a  ga 

^  ....  -  /  / 


(2.2-30) 


go 


;  i  ,  ^)  exp  (+ j  2’' Af  T  ) 


exp(-j2TAK*r  )dAfdAu. 


(2.2-31) 


i 
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Additional  transform  pairs  can  be 
follows.  With  the  use  of  Eq.  (2.2-23),  Eq.  (2 
expressed  as 


Ry(4f,Av  ;  At.Ar)  » 

-o 


or 


Rg ( Af , A^ 


and 


At, Ar) 


At, Ar) exp ( 


exp(+j2itAv^)dTdr„ 
—  — o  — o 


Rt^(T.r^  ;  At,Ar)  -  ( Af , 


or 


;  At,  Ar)  »  //  Rjj(Af,A^  ;  At,Ar)exp 
exp(-j2;rA^'r^)dAfdA^ 


■><» 


Next,  if  Eq.  (2.2-28)  is  substituted  into  Eq. 


R3(Af,A.  ;  ♦,<)  -  F^F^  Ar  ^  \  ^  '• 


obtained  as 
.2-21)  can  be 

:  At,Ar) 1 

(2.2-32) 

-j2iiAf  -  )• 

(2.2-33) 

;  At,Ar) ' 

(2.2-34) 

(+j2ffAf t)* 

(2.2-35) 

(2.2-25),  then 

it,Ar)  } 


(2.2-36) 


(2.2-41) 
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and 


RjjCAf ,A^  ;  At.Ar) 


F“^F"^{Rg(Af,Av  ;  ^,0  } 

(2.2-42) 


or 


Rjj(Af,A^  ;  At.Ar) 


Rg(Af ,Av 


*,£_)  exp  (+j2it$ At)  • 


exp(-j2itK»Ar)d»d£.  (2.2-43) 


The  various  Fourier  transform  pairs  are  summarized 
in  Fig.  6.  The  scattering  function,  or  any  of  its  Fourier 

transforms,  is  a  complete  second  order  statistical  description 
of  a  WSSUS  communication  channel. 

Another  very  important  Fourier  transform  relation¬ 
ship  can  be  obtained  by  substituting  Eq.  (2.2-35)  into  Eq. 
(2.2-28)  which  yields 


<> ,  < ) 


f,,.F^^F''b"^{R„(Af  ,Av 
At  Ar  Af  Au  H  — 


At, Ar)  } 

(2.2-44) 


or 


At,  Ar)exp(-j2:rtiAt)  • 


exp(+j2jt£- Ar)  exp(  +  j2itAfT)  • 
exp(- j2irAv"r^)dAtdArdAfdA_^ 


(2.2-45) 


and 
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since 

R„(if 


Rjj(Af,av  ;  At.ar)  =  {Rg(T,^  ;  ♦,£)) 

(2.2-46) 


Av  ;  At,Ar) 


«  oa  «  «» 

//// 


Rc  I T ,  ^ 


()!,<_)  exp  (+ j  2:t$  At)  ■ 


exp(-j2it£_'Ar)exp(-j2iiAfT)  • 
exp  ( + j  2  It  A )  d  ♦(!£ d  T  dr^ . 


(2.2-21) 
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2.2.4  Input-output  relations 

We  will  now  proceed  to  derive  an  expression  for 
the  autocorrelation  function  of  the  output  from  a  linear, 
time-varying,  space-varying,  random  filter.  The  output  auto¬ 
correlation  function  Ry (t, t* ,r,r * )  is  defined  as  follows: 

Ry<t,t‘ ,r,r*)  =  E{y(t*,r)y*(t  £')  }.  (2.2-47) 

If  Eg.  (2.1-30)  is  substituted  into  Eg.  (2.2-47),  then  the 
output  correlation  function  in  time  and  space  is  given  by 


Ry (t,t' ,r,r' ) 


•••  f  x(t-T,  r-r  ) X* (t* -T ' ,  r'-r!)* 

I  —  — o  —  -o 

exp  [+j2ir  (4t-^  '  t ')]exp  [- j  2:i  ( <_•  r  -  £'•£')]• 

Rc  (t,  T  '  ,£_,£!  ;  ,  4i ',<,<’)  dijdicdtdr* 

diji'dic'dT’dr^  .  (2.2-48) 


If  the  random  filter  is  a  WSSUS  (wide-sense  stationary  uncor¬ 
related  spreading)  communication  channel,  then  the  autocorre¬ 
lation  function  of  the  spreading  function  is  given  by  Eg. 
(2.2-16).  Substituting  Eg.  (2.2-16)  into  Eg.  (2.2-48)  yields 
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RyCt.t' .r.r'l 


x(t-T,r-r  )x*(t'-T 
—  — o 


r  *  -r 
—  -o 


)  • 


exp(+j2ii*4t)exp(-j2iT£tAr)  • 
Rg(T,r^  ;  ♦,£)di>d£dTdr^. 


(2.2-49) 


The  mean  squared  value  of  y(t,r)  is  given  by 
E{|y(t,r)|^}  »  Ry(t,t,r,r) 
where,  from  Eq.  (2.2-48),  we  have  in  general  that 


Ry(t,t,r,r) 


•  J  x(t-T,r-r^) X* (t-T ' ,r-r^) • 

exp(+j2»4*t)exp(-j2iti<_*r)  • 

R_(T,T’,r  ,£•  ;  ♦,^’,ic,<')d*dicdTdr  • 
di'dic'dr'dr*  (2.2-50) 


where  and  A£  »  In  the  case  of  a  WSSUS  cosnsu- 

nxcation  channel,  we  have  from  Eq.  (2.2-49)  that 
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Ry (t,t,r,r) 


A  oi  •  ot 

//// 


|x(t-T,  r-r^)  1 


d^djcdTdr^. 


(2. 


If  we  define  the  energy  of  the  output  signal  as 


m  m 

If 


Ey  ^  /  /  |y(t,r)  l-^dtdr,  (2. 


then  the  average  output  energy  can  be  expressed  as 


I  ^  E 
y 


«•  m 

(Ey)  -  I  j 


E{  lyu,r)  r}dtdr  (2. 


or 


m  m 

■I! 


Ey  ■  I  /  Ry(t,t,r,r)dtdr.  (2. 


If  Eg.  (2.2-51)  is  substituted  into  Eq. (2.2-54),  then 


m  m  w  m 

V'«  ■  /  /  /  /  "s' 


T,r^  '  ♦» £.) d^dKdtdr 


2-51) 


2-52) 


2-53) 


2-54) 


(2.2-55) 
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where 


«• 

j  f 


)  I  dtdr 


(2.2-56) 


is  the  energy  of  the  input  signal.  Equation  (2.2-55)  indicates 
that  the  ratio  of  the  output  (received)  average  energy  to  the 
input  (transmitted)  energy  for  a  WSSUS  communication  channel 
can  be  obtained  by  integrating  the  scattering  function  of  the 
channel.  Also  note  that  the  average  output  energy  is  not  a 
function  of  the  input  signal's  shape. 

Alternate  expressions  for  the  output  autocorrela¬ 
tion  function  can  be  obtained  from  Eq.  (2.1-25).  If  Eq. 

(2.1-25)  is  substituted  into  Eq.  (2.2-47),  then 


Ry(t,t' »£«£') 


SB  w 

I  ill 


)X*(f ' ,v') 


R„(f,f',v.v-  ;  t,f,r,r')- 
exp[+j2n  (ft-f  '  t')]exp[-j2Ti  •£' )  ] 

dfdvdf'dv’.  (2.2-57) 


If  the  random  filter  is  a  WSSUS  communication  channel,  then 
the  autocorrelation  function  of  the  transfer  function  is  given 
by  Eq.  (2.2-20).  Substituting  Eq.  (2.2-20)  into  Eq.  (2.2-57) 
yields 
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Ry(t,t' .r.r') 


X(f,v)X*lf’,v’)Rjj(af,Av 


it, ir)  • 


exp[+j2ii(ft  -  f  '  t* )  J  exp  [- j  2it 


dfdvdf ’dv* . 


(2.2-58) 


A  relationship  between  the  input  and  output  power  spectral 
densities  will  be  obtained  next. 

Let  us  assume  that  the  input  signal  x(t,r)  is  a 
zero-mean,  wide-sense  stationary  (in  time  and  space) ,  random 
process  which  is  uncorrelated  with  the  transfer  function. 
Under  these  assumptions,  the  output  autocorrelation  function 
given  by  Eq.  (2.2-57)  becomes 


,r,r') 


E{X(f,v)X*(f ’ ,y')  )• 

Rj,(f.f,v,v-  ;  t,f, £,£-)• 

exp(+j2it  (ft-f  't’)]exp  [-j2it  (v.’£-v'  •£' )  ]  • 


dfdvdfdv'  . 


(2.2-59) 


Since  x(t,£)  was  assumed  to  be  wide-sense  stationary  in  both 
time  and  space,  then  it  can  be  shown  that 


where 


E{X(f,v)X*(f’,  v’)  }  =  Sj^(f,v)6(f-f')«(y-l') 

(2.2-60) 


X(f,  v)  =  F^Fj.{x(t,£)  }, 


I 

I 


L 
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^  (2.2-61) 

is  the  power  spectral  density  of  the  input,  and 


Rjj(At,ir)  =  E{x(t,r)x*  (t*  ,r  ’ )  } 


where  it  =  t-t'  and  ir  =  .  If  Eq.  (2.2-60)  is  substituted 

into  Eq.  (2.2-59),  then 


Ry(t,f,r,r') 


S^(f ,^)Rjj(f ,f ,y^,^  ; 

exp  ( -  j  2it^-  ir )  df d^. 


t,  t',r ,  r')exp  (  + j  2:if  it)  • 
(2.2-62) 


and  if  it  is  further  assumed  that  the  transfer  function  is  wide- 
sense  stationary  in  both  time  and  space,  then  Eq.  (2.2-62)  re¬ 
duces  to 


Ry (it,ir) 


Sx(f  ,v)Rjj(f,f  ,v,  V 


it,  ir  )• 


exp(+ j2:Tf  it)  exp  (- j  2:r^-  ir)  dfd^ 


(2.2-63) 


which  implies  that  the  output  y(t,r)  is  also  wide-sense  sta¬ 
tionary  in  both  time  and  space.  If  we  define  the  output  power 
spectral  density  Sy(n,£)  as 


S.  (n,£.)  = 

y 


(2.2-64) 
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then  substituting  Eq.  t2.2-63)  inco  the  R.H.S.  of  Eq.  (2.2-64) 
yields 


Sy(n.i) 


m  «» 

I  I  (f  ,u^)  (f ,  f  ,y_,  ;  n-f,  £-y^)dfdy^ 

--  --  (2.2-65) 


where,  from  Eq.  (2.2-41), 


// 


;  ri-£,3_-v)  “II  '  Atf  ir)  exp  [-j  2- ( n-f )  it] 


exp  [  + j  2Tr  (^-y_)  •  ir  ]  ditdir . 

(2.2-66) 


Note  that  Eq.  (2.2-65)  is  in  the  form  of  a  multidimensional 
convolution  integral  analogous  to  Eq.  (2.1-46)  for  the  deter¬ 
ministic  case.  The  convolution  process  accounts  for  the  fre¬ 
quency  spreading  (frequency  in  HZ.  and  spatial  frequencies  in 
cycles/m.)  of  the  input  power  spectral  density. 


Example  2.2-1 

Using  Eqs.  (2.2-65)  and  (2.2-66),  let  us  compute  the 
output  power  spectral  density  for  the  case  when  the  transfer 
function  is  time- invariant,  space- invariant,  and  deterministic . 
First  note  that 


Rjj(f,f,v,_^  ;  It.ir)  =  E{H(f,v  ;  t,r)H*(f,v  ;  t',r')  ;  . 
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with  Eq.  (2.1-48) ] 

y  * 
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til.  COUPLING  EQUATIONS 

We  will  now  present  those  equations  which  couple  the 
transmitted  and  received  electrical  signals  to  the  transfer 
function  of  the  ocean  medium  via  the  transmit  and  receive  far- 
field  directivity  functions.  Referring  to  Fig.  2,  the  fre¬ 
quency  and  angular  spectra  of  the  input  acoustic  signal  to 
the  medium  and  the  output  electrical  signal  from  the  receive 
aperture  (array)  are  given,  respectively,  by  the  following 
equations  [2]  : 


Xj^(f ,  V) 


SB 


J X(f,a)D^(f 


v-a) da 


(3-1) 


and 


30 

y(n,r)  =  J  3)Dj^(ri,  v-e)dS 


(3-2) 


where  the  frequency  and  angular  spectrum  of  the  output  acoustic 
signal  from  the  medium  is  given  by  [see  Eq.  (2.1-46)] 


If 


n-f , 8-v) dfdv 


(3-3) 
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is  the  frequency  and  angular  spectrum  of  the  transmitted  elec¬ 
trical  signal,  and 


D,^(f,a)  =  P^{A,j,(f ,r)  )  =  J A^(f ,r)exp{+j2ira-r)dV 


(3-8) 


is  the  far-field  directivity  function  (beeun  pattern)  of  the 
complex  transmit  aperture  A^(f,r).  Also,  w.r.t.  Eg.  (3-6), 


X(f,r)  =  F^{x(t,r)  }  =  j  x(t,r)exp(-j2irft)dt 

••a* 

(3-9) 

and 

so 

A^(f ,r)  =  P~^{D^(f,a)  }  =  /  D^(f ,o)exp(-j2ita-r)da  . 

w  C0 

(3-10) 

Similarly,  w.r.t.  Sq.  (3-2),  we  have  the  following 
additional  relationships: 


Y(n,Y) 


(t,r)  } 


y  (t,r)  exp(+j2:r^-r)  dV- 


exp(-j2:trit)dt 


(3-11) 
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where  [ 2 ] 


y (t,r) 


j  yfjtn,r)Aj^(n.r)expC+j2Tint)dn, 


C3-12) 


the  frequency  and  angular  spectrum  of  the  output  acoustic  sig¬ 
nal  from  the  medium  is  given  by 


yj^(t,r)  exp(+j2n^- 1)  dV* 
exp(-j2unt)dt  (3-13) 


where  [see  Eq.{2.1-25)] 


QD  40 

yj^(t,r)  =  j  j  Xj^(£,^)Hj^(f ,y^  ;  t,r) exp (+ j2nf t)  • 


exp(-j2-«r\^*r)dfdy_  , 


(3-14) 


and 


Djj(n,£)  =  F^{Aj^(n,£)}  =  J  Aj^(n.r)exp(+j2TT£'r) 


dV 


(3-15) 


is  the  far-£ield  directivity  function  (beam  pattern)  of  the 
complex  receive  aperture  Ajj(n»r).  Finally,  w.r.t.  Eq.  (3-12), 
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a» 

“  J  yM(^t,r)exp(-j2itn  t)dt 


(3-16) 


and 


Aj^(n.r)  =  F“^{Dj^(n  ,8.)}  =  J  Dj^  (n  ,8.)  exp  (- j  2ii  3.‘ r)  d^. 

■■CO 

(3-17) 


Substituting  Eg.  (3-14)  into  Eq.  (3-16)  yields 


40  CO 

Yjj(n,r)  j"  '  '^"f/£.)8xp(-j2nv*£)dfd^ 


_dO  .00 


(3-18) 


where 


Bj^(f,y_  ;  ♦,£)  =  J  Hj^(f,y.  ;  t,r)exp(-j2Ti((it)dt 
*00 

(3-19) 

which  is  not  the  same  as  the  bi-frequency  function  [see  Eq. 

(3-4) ]  . 

The  components  of  the  vectors  a,  £,  and  £  are  spatial 
frequencies  with  units  of  cycles/m.  and  r  =  (x,y,z). 

A  moments  reflection  leads  one  to  the  conclusion  that 
Eq.  (3-12)  is  well  suited  for  space-time  signal  processing 
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applications.  By  inspecting  Eqs.  C3-18) ,  C3-19) ,  and  (3-1), 
it  can  be  seen  that  the  output  signal  given  by  Eq.  (3-12)  can 
be  expressed  as  a  function  of  tine  ^md  space  in  terms  of  the 
transmitted  electrical  signal,  the  transmit  aperture,  the 
transfer  function  of  the  ocean  medium,  and  the  receive  aperture. 
Equation  (3-2) ,  however,  represents  the  theoretical  output 
"spectrum",  albeit  dependent  upon  the  same  system  functions. 

If  the  ocean  medium  is  modelled  as  a  random  filter, 
then  y(t,r)  given  by  Eq.  (3-12)  is  a  random  process  with  auto¬ 
correlation  function 

Ry(t,t',r,r')  =  E{y(t,r)y*(t’,r’) ).  (3-20) 

It  is  probably  obvious  by  now  that  the  weak  link  in  the 
coupling  equations  is  the  ocean  medium  transfer  function 
Hjj(f,^  ;  t,r)  and  its  corresponding  autocorrelation  function 
[see  Sq.  (2. 2-2) ] 

R„  (f,f,v,.'  ;  t,f,r,r')  »E{H„(f,v  ,  t,£)H*(f,.'  ;  f,r')}. 

M 

(3-21) 

If  one  is  not  able  to  specify  a  realistic  functional  form  for 
Hjj(f,^  ;  t,r)  ,  then  the  equations  presented  in  this  section 
represent  an  exercise  in  the  methods  of  linear  system  theory. 
This  leads  us  to  the  main  purpose  of  this  paper  which  is  to 
demonstrate  the  derivation  of  an  ocean  transfer  function  based 
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upon  the  W.K.B.  approximation.  The  derivation  is  discussed  in 
Section  4.1. 

Finally,  note  that  the  autocorrelation  function  given 
by  Eq.  (3-21)  can  be  referred  to  as  a  generalized  coherence 
function  since  it  is  a  generalization  of  the  two-frequency  cor¬ 
relation  function  or  two-frequency  mutual  coherence  function 
based  upon  linear,  time-variant  filter  theory  as  discussed  by 
Ishimaru  [51],  for  example.  Equation  (3-21)  provides  informa¬ 
tion  concerning  the  2uiiount  of  spreading  in  time-delay,  space, 
frequency,  and  spatial  frequencies  that  a  transmitted  signal 
will  be  expected  to  experience  as  it  propagates  in  a  random 
medium.  The  result  of  all  this  spreading  is,  of  course,  dis¬ 
tortion  in  pulse  shape. 
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IV.  ANALYSIS 

4.1  Transfer  Function 

We  will  now  proceed  to  derive  a  transfer  function 
which  models  the  bistatic  communication  channel  geometry  shown 
in  Fig.  1.  The  communication  channel  is  regarded  to  be  the 
ocean  volume  between  the  apertures  so  that  surface  and  bottom 
scattering  effects  are  not  included.  Both  apertures  are  sta¬ 
tionary  (not  in  motion) ,  and  it  is  assumed  that  no  discrete 
point  scatterers  (such  as  bubbles,  fish,  etc.)  are  in  the 
volume  between  the  apertures.  No  motion  implies  that  the  re¬ 
sulting  transfer  function  will  be  time-invariant. 

The  propagation  of  small  amplitude  acoustic  signals 
in  the  ocean  from  the  transmit  aperture  to  the  receive  aperture 
can  be  described  by  the  following  linear,  inhomogeneous,  scalar 
wave  equation: 


7‘‘v»(t,r)  - 


c^(r) 


3t 


J  f(t,r)  a  Xjj(t,r) 


(4.1-1) 


where  F(t,r)  is  the  velocity  potential  at  time  t  and  position 
r  -  (x,y,z),  Xj^(t,r)  is  the  source  distribution  [see  Fig.  2 
and  Eq.  (3-6)],  and  c(r)  is  the  speed  of  sound  in  the  ocean. 
Since  the  coupling  equations  discussed  in  Section  III  already 
allow  for  an  arbitrary  x„(t,r)  with  corresponding  frequency 
and  angular  spectrum  Xj^(£,y^) ,  we  need  only  find  the  solution 
to  the  following  Helmholtz  wave  equation: 


62 


7^vCrl  +  n^(r)¥>(r)  =  0  C4.1-2) 

where 

=  2iif/e  =  2n/X  (4.1-3) 

o  o  o 

is  the  constant,  reference  wavenumber, 

n(r)  =  Cg/cCr)  (4.1-4) 

is  the  random  index  of  refraction, 

(4.1-5) 

o  — o  o 

is  the  constant,  reference  speed  of  sound  at  the  source  posi¬ 
tion  r^  =  ^^o'^o'^o’ ' 

<P{t,r)  =  ^(r)exp(+j2iTft)  (4.1-6) 

is  the  time-harmonic  solution  of  Eq.  (4.1-1)  when  x„(t,r)  is 
set  equal  to  zero, and  where  i?(r)  is  the  solution  of  Eq.  (4.1-2) 
Note  that  the  wavenumber 

k(r)  =  2iff/c(r)  (4.1-7) 


can  be  expressed  as 


k(r)  =  • 


(4.1-8) 


L 
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Therefore,  k(r  )  =  since  n(r  )  =  1.  The  index  of  refraction 
— o  o  —o 

is  commonly  written  as  [40,  42-45] 

n(r)  =  njj(r)  +  n^^Cr)  (4.1-9) 

or 

n(r)  =  njj(r)  +  a(r)n^j^(r)  (4.1-10) 

where  nQ(r)  is  the  deterministic  component  and  is  usually 
close  to  unity  in  value,  is  the  random,  zero  mean  compo¬ 

nent,  air)  is  the  standard  deviation  of  njj(r)  ,  and 


njjj^(r)  =  nj^(£)/a(r)  (4.1-11) 

is  the  normalized  random  component  with  zero  mean  and  variance 
equal  to  unity.  We  shall  work  with  Eq.  (4.1-10)  in  this  paper. 
Note  that  the  average  value  of  n(r)  is  equal  to  njj(r). 

Let  us  assume  that  the  speed  of  sound  is  only  a 
function  of  the  depth  y,  i.e.,  c(r)  =  c(y),  so  that  Eq.  (4.1-2) 
reduces  to 

V^'^(x,y,z;  +  k^  n^(y)<^(x,y,z)  =  0  (4.1-12) 

where,  from  Eq.  (4.1-10), 

n(y)  = 


np(y)  +  o(y)n^j^(y)  . 


(4.1-13) 
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By  using  the  method  of  separation  of  variables  and  the  W.K.3. 
approximation  [46,47],  an  approximate  solution  of  Eq.  14.1-12) 
is  given  by 


<P(x,y,z)  a  exp[-jkjj(x-XQ)]  (y)  exp[+ je  y  (y)  ]  exp[  - jk^  ( z-z^)  ] 

(4.1-14) 


where 


aY(y) 


[k^Cy)] 


-1/2 


(4.1-15) 


9y<y) 


exp 


ky  {;)di; 


(4.1-16) 


and 


2.2, 


K^iy)  =  [k^n  (y)  -  "z 


k^-k^] 


1/2 


(4.1-17) 


where  k^,  k^(y).  and  k^  are  the  components  of  the  propagation 
vector 


k(y)  =  kjjX  +  k^Cy'y  +  k^z.  (4.1-18) 

Note  that  k..  and  k_  are  constants  while  k^Cy)  is  a  function  of 
X  L  ^  ^ 

the  depth  y.  In  addition,  Eqs.  (4.1-14)  thru  (4.1-16)  allow 

for  a  general  source  location  with  position  vector  =  (x^,y^,z^). 

The  W.K.B.  approximation  given  by  Eqs.  (4.1-15)  and 
(4.1-16)  is  a  valid  solution  for  the  depth  dependence  provided 
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that  [46,47,50]  1)  the  transmitted  frequencies  f  are  high,  2) 
the  sound  speed  profile  c(y)  is  slowly  changing,  and  3)  the 
depth  interval  from  y^  to  y  does  not  include  any  turning 
points.  A  turning  point  exists  at  y  =  y^  if  kyCyip)  *  0. 

Since  ^he  reference  propagation  vector  can  be 

expressed  as 


=  k^x  +  Ic^y  +  k^z  ,  (4.1-19) 

then 

•^o  =  IlSol^  =  ^  ^  ’'Z  (4.1-20) 

where 


k 


X 


k  u 
o  o 


(4.1-21) 


and 


where 


u 

o 


=  sine 


k  V 

o  o 


k  w 
o  o 


^cos^o 


(4.1-22) 


(4.1-23) 


(4.1-24) 


v^  =  sineQSi.iii/^  =  cosBg  (4.1-25) 

and 

w^  =  cose^j  (4.1-26) 

are  the  direction  cosines  w.r.t.  the  positive  X,y,  and  Z  axes, 
respectively,  and  (Sgf'l'g)  are  the  vertical  and  azimuthal 
spherical  angles  measured  w.r.t.  the  positive  Z  and  X  axes, 
respectively,  representing  the  initial  directions  of  wave  pro- 


i 
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pagation  (see  Fig.  7).  Prom  Eq.  (4.1-20)  we  can  write  that 


(4.1-27) 


and  substituting  Eqs.  (4.1-27)  and  (4.1-22)  into  Eq.  (4.1-17) 
yields 

(  .1/2 


=  ky 


1  ^  [nNy)-ll 


(4.1-28) 


Note,  that  since  n(yg)  =  1.  Therefore,  the  square 

root  expression  in  Eq.  (4.1-28)  is  responsible  for  changing 
the  initial  direction  of  propagation  )tY.  If 

i [n^(y)  -  1]/Vg|  <  1  ,  (4.1-29) 

then  the  square  root  expression  in  Eq.  (4.1-28)  can  be  approx¬ 
imated  by  the  first  two  terms  in  a  binomial  expansion  yielding 

k.,(y)  w  ky  +  k^[n^(y)  -  1]/ [2^^)  .  (4.1-30) 

It  will  be  shown  later  that  the  binomial  expansion  criterion 
given  by  Sq.  (4.1-29)  can  be  related  to  the  critical  angle  of 
incidence,  and  hence,  total  reflection.  The  assumption  con¬ 
cerning  the  absence  of  turning  points  will  also  be  discussed 
later,  but  first  let  us  return  to  the  derivation  of  the  trans- 


er  function. 
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Recall  that  we  designated  the  vector  ^  to  represent 


the  transmitted 

spatial 

frequencies  as  in  X„(f,v) 

(see  Fig.  2 

If  we  let 

(4.1-31) 

then 

= 

(4.1-32) 

fy 

=  kY/(2ir)  =  v^/\^ 

(4.1-33) 

and 

^2 

=  k2/(2Tt)  = 

(4.1-34) 

are  the  spatial  frequencies  in  the  X,Y,  and  Z  directions,  re¬ 
spectively,  where  c^  =  Therefore,  if  Eqs.  (4.1-30) 

thru  (4.1-34)  are  substituted  into  Eqs.  (4.1-14)  thru  (4.1-16), 
then  Eq.  (4.1-6)  becomes 


e(t,x,y,z)  « 

j2,rf^  +  [n^(y)  -  1]  /(ATtf^)! 

exp  (+ j2itf  t)  exp  [- j2ii^*  (r-r^)  ] .  (4.1-35) 

Since  the  input  to  the  communication  channel  is  the  time-har¬ 
monic  plane  wave 

exp  (+j2iTft)  exp  [- jZuv"  (£-£q)  ]  , 
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and  since  Eq.  (4.1-3S)  represents  the  output  from  the  channel 
at  any  tine  t  and  position  r,  then  [see  Eq.  (2.1-21)] 

(4.1-36) 


where  the  random,  time-invariant,  space-variant  transfer  func¬ 
tion  of  the  ocean  medium  is  given  by 


where 


;  y)  =  '■  y>e5cp[+j0j^(f,fy  ;  y)] 

(4.1-37) 


Ajj(f,fY  ;  y)  =  {2itfY  +  lc^[n^(y)-l]/(4irfy) 

(4.1-38) 


3j^(f,fy  ;  y)  =  -[k^/(4„fY)  ]  /  [n^(c)-l]dc 


'o''  “'■Y' 


(4.1-39) 


■"o  “ 


(4.1-3) 


and  from  Eq.  (4.1-13), 

n^(y)  =  np(y)  +  2np  (y)  a  (y)  (y)  +  (y)  n^P  (y)  . 

(4.1-40) 


If  the  medium  is  "weakly  irregular"  or  "weakly  scattering" , 
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i.e.,  if  a (y)  is  very  small  compared  to  unity  so  that 


<<  1 


(4.1-41) 


where  s^  is  the  scale  size  of  the  irregularities  (i.e.,  the 

average  distance  over  which  the  refractive  index  fluctuations 

2 

remain  correlated  [48]),  then  terms  involving  a  (y)  can  be  ne¬ 
glected  [42]  and  Eq.  (4.1-40)  reduces  to 


n^(y)  »  b^(y)  +  2np(y)o  (y)nj,jj(y) .  (4.1-42) 


Therefore,  by  using  Eq.  (4.1-29)  to  simplify  Eq.  (4.1-38)  and 

* 

upon  substituting  Eq.  (4.1-42)  into  Eq.  (4.1-39),  one  obtains 
the  following: 


Aj^(f,fY  ;  y)  “  {2Trfy) 


-1/2 


(4.1-43) 


and 


(4.1-44) 

where  y 


6„i,(f,fY  ;  y)  -  -[k‘/ 


(^itfy)]  J 


[n^(c)-l]dc  (4.1-45) 


is  the  deterministic  or  average  component  of  the  phase  function. 


and 
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y 

^o  (4.1-46) 

is  the  random  coa^nent. 

The  medium  transfer  function  given  by  Eq,  (4.1-37) 
with  amplitude  and  phase  functions  and  6^^  specified  by  Eq. 
(4.1-43)  and  Eqs.  (4.1-44)  thru  (4.1-46),  respectively,  in¬ 
dicate  that  for  a  weaJtly  irregular  medium,  the  major  effect  of 
the  medium  is  to  angle  modulate  the  transmitted  field.  Further¬ 
more,  If  njj(y)  «S  1,  then  aj^(f,fy  ;  y)  a  0  [see  Eq.  (4.1-45)] 
and,  as  a  result,  the  angle  modulation  is  due  strictly  to  the 
random  fluctuations  o  (y)nj,j^(y)  of  the  index  of  refraction  [see 
Eq.  (4.1-46)].  Note  that  the  angle  modulation  process  is  often 
referred  to  as  "scattering"  [48].  Also  note  that  the  transfer 
function  derived  in  this  section  can  be  written  as  the  product 
of  two  functions,  one  deterministic  and  the  other  random,  i.e., 

Hm  “  [exp(+j0jy^)] 

which  agree!'  with  the  assumed  transfer  function  expressions  of 
Laval  and  Labasque  [36]  and  with  the  general  practice  of  re¬ 
presenting  a  field  in  a  random  medium  as  the  product  of  a  de¬ 
terministic  and  a  random  function  [40]. 

Let  us  now  discuss  the  two  important  assumptions 
responsible  for  the  derivation  of  the  transfer  function,  namely. 
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the  aforementioned  binomial  expansion  criterion  given  by  Eq. 
(4.1-29)  emd  the  assumption  regarding  the  edisence  of  turning 
points. 

Since 

n(y)  »  c^/c(y)»  (4.1-47) 


substituting  Eqs.  (4.1-25)  auid  (4.1-47)  into  Eq.  (4.1-29) 
yields 


8  <  cos 

o 


-1 


(Ic2(y)  - 


c(y> 


(4.1-48) 


where  8^  is  the  angle  of  incidence  of  the  reference  propagation 
vector  (see  Fig.  7) .  Now  recall  that  the  critical  angle  of 
incidence  8^  associated  with  a  time-harmonic  plane  wave  inci¬ 
dent  upon  a  plane  boundary  between  two  fluid  media  is  given  by 
[49]  (see  Fig.  8) 


sinS^j  “  °l/®2  '  *^2  ^  *^1  (4.1-49) 


or,  equivalently. 


8  *  cos 

c 


Cj  >  (4.1-50) 


where  Cj  must  be  greater  than  Cj^  for  8^  to  exist.  When  the 

angle  of  incidence  8  >  8  ,  there  is  total  reflection  and  hence, 

c  — — .  ■ 


74 


no  transmission  of  energy  into  medium  II.  The  speed  of  sound 
c(y)  generally  increases  with  depth,  for  example,  in  the  deep 
ocean.  Also,  cCyl  increases  above  and  below  the  SOFAR  channel 
axis.  Therefore,  when  c(y)  >  c^,  the  absolute  value  sign  in 
Eg.  (4.1-48)  can  be  removed,  and  by  comparing  Eg.  (4.1-50) 
with  the  R.H.S.  of  Eg.  (4.1-48),  the  binomial  expansion  crite¬ 
rion  indicates  that  the  initial  angles  of  transmission  must 
be  less  than  the  "critical  angle"  in  order  to  avoid  total  re¬ 
flection,  and  thus,  passing  thru  a  turning  point  [47,50]. 

This  is  consistent  with  the  fact  that  the  W.R.B.  approximation 
is  invalid  at  a  turning  point  [47,50]. 
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4.2  Output  Electrical  Signal 

Now  that  we  have  derived  an  ocean  medium  transfer 
function,  let  us  demonstrate  the  use  of  the  coupling  equations 
presented  in  Section  III  by  calculating  the  output  electrical 
signal  y(t,r)  from  the  receive  aperture  (array)  as  given  by 
Eq.  (3-12) .  As  can  be  seen  by  inspecting  the  integrand  of 
Eq.  (3-12),  we  need  expressions  for  both  of  the  kernels  Yjj(n,r) 
and  Aj^(n,r).  Let  us  compute  y(j(n,r)  first  [see  Eq.  (3-18)]. 

Assume  that  the  transmit  aperture  depicted  in  Fig.l 
is  a  planar  array  of  MxN  (odd)  complex  weighted  point  sources, 
centered  at  parallel  to  the  XY  plane.  In  addi¬ 

tion,  assume  that  the  complex  weights  are  separed^le.  Since  in 
most  practical  situations  an  identical  input  electrical  signal 
is  applied  to  all  elements  in  the  transmit  array  before  the 
complex  weights,  i.e.,  since 

x(t,r)  H  x(t) , 

then  Eq.  (3-7)  reduces  to 

X(f,a)  =»  X(f)  5(a)  (4.2-1) 

since  F^{1}  «  5  (a^)  .  Substituting  Eq.  (4.2-1)  into  Eq.  (3-1) 
yields 

Xj^(f,v)  -  X(f)D^(f,v)  (4.2-2) 
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where 

(M-l)/2  (N-l)/2 

Cj|jd^exp(+j2irfj^djj)  • 

m— CM-l)/2  n»-(N-l)/2 

exp(+j2iTfYnd^)exp(+j2irfjjX^)  • 
exp(+j2TrfYyQ)exp(+j2*f2ZQ) 
(4,2-3) 

is  the  far-field  beam  pattern  of  the  transmit  array,  c_  and  d 

m  n 
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m 

Yj^(n  ,r)  “  X(n  )  /  )Hjj(n,v  ;  r)exp(-j2i'^*r)d^, 

■*  (4.2-6) 

and  upon  substituting  Eg.  (4.2-6)  into  Eg.  (3-12),  we  obtain 


y(t,r)  =  j  X(f) 

Aj^(f  ,r)exp(+32iift)df  (4.2-7) 

where  n  was  replaced  by  £  since  there  is  no  frequency  (HZ) 
spreading  in  this  example  due  to  tlie  time- invariant  property 
of  the  transfer  function.  We  need  to  specify  the  complex 
aperture  function  Ap(f,r)  next. 

Assvime  that  the  receive  aperture  depicted  in  Fig.  1 
is  a  planar  array  of  M'x  N’  (odd)  complex  weighted  point 
sources,  centered  at  (Xj^,y^,z^)  and  parallel  to  the  XY  plane. 

In  addition,  assvune  that  the  complex  weights  are  separable. 
Therefore,  the  receive  aperture  function  is  given  by 


J"  '  r)exp(-j2iT^*r)d^’ 


Aj^(f  ,r) 


(M'-l)/2  (N’-l)/2 

£  I]  +  • 

i— (M’-l)/2  g— (N’-l)/2 

«(y“tyR  +  gd^] )  5  (2-Zj^) 


(4.2-8) 
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where  c|  and  are  complex  weights  and  d^  and  d^  are  the  inter¬ 
element  spacings  in  the  X  and  Y  directions,  respectively. 

Upon  substituting  Eqs.  (4.1-31),  (4.1-36),  (4.2-3), 

and  (4.2-8)  into  Sq,  (4.2-7)  and  recalling  that  r  =  (x,y,z), 
one  obtains 

(M'-l)/2  (N'-l)/2 

y(t,x,y,z)  =2]  cM^  f  X(f)- 

i=-(M'-l)/2  q=-(N‘-l)/2 

(M-l)/2 

Z  ^m  f 

m»-(M-l)/2  -» 

(N-l)/2 

“^n  (  •  V  qdy)exp(-j2^fylY^^)dfY- 

n=-(N-l)/2  -<• 

00 

J  exp(-j2'!rf2iZ)df2exp(  +  j2:rft)df • 

mOO 

6(x-[Xj^+  idj^] )  «  (Y-[yjj+  qd^] )  6  (z-Zj^)  (4.2-9) 

where 

^^im  =  ^  (4.2-10) 

■  "‘^Y^  (4.2-11) 

and 


AZ  =  Zj^  -  z^. 


(4.2-12) 
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If  we  now  change  varicUsles  from  spatial  frequencies  to  direc¬ 
tion  cosines  by  substituting  Egs.  (4.1-32)  thru  (4.1-34)  into 
Eq.  (4.2-9),  and  treat  the  frequency  variable  f  as  a  constant 
with  respect  to  the  spatial  frequency  integrations,  then  Eq. 
(4.2-9)  becomes 


(tr-l)/2  (N'-l)/2 


y(t,x,y,2)  =  {2/c^)  ^ 


i=-(M'-l)/2  q— (tJ’-l)/2 
(M-l)/2 


E  '1^4  / 


f  X(f)  • 


E 

m*-  m~l)/2 


c„sinc(2fAX.  /c J 
tn  im  o 


(N-l)/2  +1 

n=-(N-l)/2  a^ 
exp(-j2^[fv^/c^]aYg^)dv^j. 

{b„exp(-jiib  faz/c  )sinc(b  faz/c  )}• 

H  H  o  C5  o 

exp(+j2iTft)df  6(x-  [Xjj+  idj^])S(y-  [yj^+  qd^] )  5  ( z-z^^) 


(4.2-13) 


since  in  our  exzunple  problem 
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-1  <  <  +1 
o 


(4,2-14) 


and 


where 


a  <  V_  <  +1 

C4.2-15) 

q  o 

0  <:  w  <  b_ 

(4.2-16) 

-  -  O  q 

^<5= 

(4.2-17) 

bq  =  (1  -  ^ 

(4.2-18) 

,  .  sin(irx) 

3inc(x)  = 

(4.2-19) 

and  ( • )  was  replaced  by  nQ(*)  in  Eqs.  (4.2-17)  and  (4.2-18). 
The  transfer  function  is  given  by  Eqs.  (4.1-37)  and  (4.1-43) 
thru  (4.1-46)  with  the  exception  that  the  lower  limit  of  in¬ 
tegration  y^  in  Eqs.  (4.1-45)  and  (4.1-46)  must  be  replaced  by 
V  +  nd„  which  is  the  Y  coordinate  of  a  point  source  in  the 

“  O  Y 

transmit  planar  array.  Note  that  if  njj(*>  <  1,  then 


Equation  (4.2-13)  is  the  desired  result,  i.e.,  it  represents 
the  random  output  electrical  signal  from  each  element  in  the 
receive  array  in  terms  of  the  transmitted  electrical  signal, 
the  transmit  and  receive  arrays,  and  the  random  transfer  func¬ 
tion  of  the  ocean  medium.  If  y(t,x,y,2)  given  by  Eq.  (4.2-13) 
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is  complex,  then  simply  take  the  real  part  to  obtain  the  real 
output  electrical  signal. 
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4 . 3  Coherence  Function 

Upon  inspecting  Eq.  (4.2-13),  it  can  be  seen  that 
the  autocorrelation  function  of  the  output  electrical  signal 
will  depend  upon  the  autocorrelation  function  of  the  transfer 
function 


M 


E{Hj^(f,f.^  ;  y)H*(f’,f^  ;  y')  } 

(4.3-1) 


which  is  also  known  as  the  coherence  function.  Substituting 
Eqs.  (4.1-37),  (4.1-39),  and  (4.1-43)  into  Eq.  (4.3-1)  yields 


Rg  (f,f',fy,f^  ;  y,y’)  *  [ait  (f^f^)  ^  sjaxp  (+ j  [KI  (y) +K  ’  I  (y  ‘ )  J ) 

M  I 


where 


K  =  -  k^/(4iif.^) 


(4.3-2) 


(4.3-3) 


K’  =  +(k^)  V(4iif^) 


(4.3-4) 


[n^(5)  -  l]di; 


[n^(0  -  l]dc 


(4.3-5) 


(4.3-6) 
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k  «  2irf/c  (4,3-7) 

o  o 

k*  -  2itf/c^  (4,3-8) 

o  o 

2 

and  n  (■)  is  given  by  Eg,  (4,1-42)  where  the  expectation  appear¬ 
ing  in  Eg,  (4,3-2)  is  the  characteristic  function  of  the  random 
guantity  [KI(y)  +  K'Ky*)]-  If  it  is  assumed  that  I  (y)  is  a 
real  Gaussian  random  process,  which  implies  that  the  index  of 
refraction  is  a  real  Gaussian  random  process,  then  Eg,  (4,3-2) 
can  be  written  as 


Rj,  (f,f',fy,f^  ;  y,y')  *  [2ii(fyf^)^/2]  exp(+j  [e„jj(f  ,f^  ,  y)- 

M 

exp(-E{e^{f,fY  ;  y)  }/2)  • 
exp(+E{0.^(f  ,fY;  y)  ej^j^(f  ’  ,f^  ;  y')))- 
exp(-E{e^(f  ,f^  ;  y’)//2)  (4,3-9) 


where 

E{9jj^(f,fY  ;  y)} 


[k2/{2irf^)]^ 


y  y 


njj(c)np(c')a(c)a(;')  • 

R_  (c,c')dcdc' 

"nr 


(4,3-10) 
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y’  y’ 


E{e^(f',f^  ;  y')}  =  [(k^)  V(2irf^)]^  j  j  ( c)  ( ?  ' )  a  (? )  a  (c - 


y  y 
•'o  ■'o 


R  (c,C’)dcdc’ 
"nr 


(4.3-11) 


y  y’ 

^  J"  njj(c)nQ(i;')o(c)a(c') 


R  (5,;')dcd;' 

"nr 


(4.3-12) 


and 


Rn^^(y.y)  -  E{n^^(y)n^^(y<)}.  (4.3-13) 


If  it  is  further  assigned  that  the  deterministic  component  n^  (y) 
of  the  index  of  refraction  is  equal  to  unity  [40,  42-45],  and 
that  the  random  component  nj^(y) »  and  hence,  the  normalized 
random  component  Ojjjjfy)  is  wide-sense  stationary,  i.e.. 


R_  (y.y')  =  R„  (4y) 
"nr  "nr 


where  Ay  ■  y-y*,  then  [see  Eq.  (4.1-45)] 


(4.3-14) 


(4.3-15) 
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and 


(4.3-16) 


and  Eqs.  (4.3-10)  thru  (4.3-12)  become,  respectively, 


1  Tk^a  f 

(y-y-) 

r  r  1,1  1 

[2ff7j 

/  [' 

(t)dc 


-(y-Yo) 


(4.3-17) 


2  (y’-y„) 


,  y)[  -  ,y-y^,[^]  j  [i  - 


-(y'-y^) 


(4.3-18) 


and 


®  i®MB  ^Y' 


(kok^o)'- 


(y-yo> 


(y-y') 


/..  ' 


•y-Jo' 


(i:)dc  +  (y*-y  ) 

"nr  ° 


(y-y') 


/ 


R  (i)dc  + 

"nr 


-(y’-Yo) 


U 

/ 


(y-y^) 


-(y'-Yo) 


?R„  (OdC  -  /  SR  (S)dS 

"nr  J  "nr 

(y-/') 


(4.3-19) 
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where  a  is  the  constant  standard  deviation  of  the  Gaussian, 
zero  mean,  wide-sense  stationary,  random  component  of 
the  index  of  refraction. 
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V.  SUMtlARY  and  COMCLUSIONS 

A  consistent  notation,  fundamental  input-output  relations, 
and  various  time-space  transformations  for  both  deterministic 
and  random  linear,  time-variant,  space-variant  filters  have 
been  established.  The  notation  is  consistent  in  the  sense  that 
all  of  the  various  input-output  relations  which  are  based  upon 
the  general  theory  will  reduce  to  the  classical  relations  of 
linear,  time-invariant  filter  theory.  These  results  should  be 
of  interest  to  persons  involved  in  the  general  area  of  linear 
systems  theory,  and  not  only  to  those  involved  in  underwater 
acoustics . 

With  the  use  of  the  method  of  separation  of  variables 
and  the  W.K.B.  approximation,  a  mathematical  expression  of  a 
time-invariant,  space-variant,  random  transfer  function  of  the 
ocean  medium  was  derived.  The  transfer  function  was  time-in¬ 
variant  instead  of  time-variant  because  motion  was  not  consi¬ 
dered  in  the  present  derivation.  The  transfer  function  modelled 
the  ocean  volume  between  transmit  and  receive  apertures  (arrays) . 
The  ocean  volume  was  characterized  by  a  random  index  of  refrac¬ 
tion  (sound  speed  profile)  which  was  a  function  of  depth.  The 
index  of  refraction  was  decomposed  into  deterministic  and  random 
components. 

In  addition  to  the  transfer  function  derivation,  two 
ex^unple  calculations  were  made.  The  first  example  demonstrated 
the  use  of  the  coupling  equations  and  involved  the  derivation 
of  a  mathematical  expression  for  the  random  output  electrical 
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signal  at  each  element  in  a  receive  planar  array  of  complex 
weighted  point  sources.  The  output  signals  were  expressed 
in  terms  of  the  frequency  spectrum  of  the  transmitted  electrical 
signal,  the  transmit  and  receive  arrays,  and  the  previously 
derived  transfer  function  of  the  ocean  medium.  The  first  exam¬ 
ple  demonstrated  that  an  output  electrical  signal  could  be 
derived  in  a  logical  and  straightforward  fashion.  The  second 
example  involved  the  derivation  of  the  coherence  function,  i.e., 
the  autocorrelation  function  of  the  transfer  function.  In 
order  to  obtain  somewhat  simplified  results,  it  was  necessary 
to  assume  that  the  random  component  of  the  index  of  refraction 
was  Gaussian  and  wide-sense  stationary. 
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